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Assessment  of  the  AN/FPS  118  Ionospheric  Model 

and  Proposed  Improvements 


1.  INTRODUCTION 


For  the  frequency  management  of  the  Over-the  Horizon  Backscatter  (OTH-B)  radar  and 
for  the  target  coordinate  registration  and  propagation  assessments.  Millman  et  al.  ^ 
developed  the  AN/FPS  118  -  Analytical  Model  Specification  henceforth  referred  to  as  the 
OTH  ionospheric  model.  The  model  is  used  to  generate  ionospheric  profiles  and  parameters 
at  the  grid  points  in  100  nautical-mile  increments  along  the  radial  from  the  radar,  covering 
the  OTH  coverage  area.  The  coordinate  registration  scheme  computes  the  virtual  height 
traces  (synthetic  ionograms)  at  the  midpoints  of  the  selected  ground  ranges  from  the  OTH 
ionospheric  model.  Superposing  the  transmission  curve  on  the  synthetic  ionogram  traces 
provides  the  reflection  height  tables  for  the  respective  frequencies.  The  derived  virtual 
height  traces  are  also  used  to  predict  operational  frequencies,  especially  at  times  and  in 
sectors  where  the  backscatter  ionograms  do  not  provide  clear  signatures  of  the  ground 
scatter. 


(Received  for  Publication  1  March  1994) 

1  Millman,  G.H..  Bowser,  C.A.  and  Swanson.  R.W.  (1988)  An  ionospheric  model  for  HF  sky 
wave  backscatter  radar,  NATO-AGARD  Symposium  on  Ionospheric  Structure  and  Variability 
on  A  Global  Scale  and  Interactions  with  Atmosphere,  Magnetosphere,  Munich,  Federal 
Republic  of  Germany.  AGARD  Conf.  Proc.  441:pp.  43-1  to  43-15. 
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The  OTH  ionospheric  model  is  an  updated  version  of  the  AFGWC  polar  model!  1982) 
developed  by  the  U.  S.  Air  Force  Global  Weather  Central  (AFGWC)  .2  The  model  uses  the 
ITS-78  ionospheric  model  of  Jones  et  al .3  as  the  starting  point.  The  basic  model  is  driven  by 
an  'effective  sunspot  number  SSN^jy’  computed  by  AFGWC  from  the  foF2  data  for  the  past  five 
days  from  a  network  of  50  ionospheric  stations.  Assuming  persistence  as  the  best 
prediction  available,  the  least  squared  deviation  fit  of  the  ionospheric  model  to  the  past  five 
days'  foF2  averages  for  each  hour,  produces  an  effective  sunspot  number  SSN^d- .  which  is 
used  to  predict  the  ionospheric  state  for  the  current  day.  Using  the  effective  auroral  oval 
parameters  Qg  and  Kp,  the  AWS  superposes  the  auroral  E  and  F  layer  enhancement  and  the 
midlatitude  F  layer  trough  electron  density  depletion  on  the  lTS-78  model.  The  ionospheric 
parameters  obtained  by  the  automated  Air  Weather  Service  (AWS)  Digital  Ionospheric 
Sounding  System  (DISS)  net  are  used  for  real-time  updating  of  the  model  at  the  radar  sites. 

The  OTH  ionospheric  model  was  tested  and  evaluated  at  the  Phillips  Laboratoiy 
(formerly  called  the  Geophysics  Laboratoiy)  after  an  in-house  capability  to  run  the  model 
emd  to  display  various  model  outputs  was  developed.  The  results  of  these  tests  along  with 
necessary  and/or  beneficial  modifications  are  presented  below. 


2.  THE  OTH  IONOSPHERIC  MODEL 


The  OTH  ionospheric  model  Implemented  in  the  East  Coast  Radar  System  (ECRS)  has 
been  documented  in  an  unpublished  paper  (private  communication.  Carl  Bowser.  GE.  1988). 
Mlllman.  Bowser  and  Swanson  (1988)^  have  published  the  model  without  the  updating 
scheme  as  a  NATO-AGAFID  symposium  report^  (Millman  et  al.  1988).  Both  reports  and  a 
computer  program  of  the  model  are  used  for  testing  the  concepts,  consistency  and  the 
working  accuracy  of  the  OTH  ionospheric  model.  A  copy  of  the  OTH  ionospheric  model 
computer  program  was  provided  to  PL  by  GE  and  has  been  installed  at  the  PL  computer 
facility.  At  the  radar  site,  the  computer  program  runs  in  a  real  time  operational  mode.  For 
analysis,  the  program  was  modified  at  PL  to  run  in  a  batch  mode  for  the  purpose  of  testing 
the  model  at  PL.  It  is  now.  for  example,  possible  to  provide  diurnal  variations  of  a  selected 
set  of  ionospheric  parameters  at  points  of  interest,  to  permit  comparison  with  test  data. 

Modifications  found  necessary  were  made  in  one  copy  of  the  model,  while  an  unmodified 
version  is  being  maintained  as  a  reference  . 


2  Air  Force  Global  Weather  Central  (1982)  AFGWC  Polar  Ionospheric  Model  Air  Force 
Weather  Central,  Program  Listing,  Feb.  1982. 

3  Jones,  W.B..  Graham.  R.P.  and  Leftin,  M.  (1969)  Advances  in  Ionospheric  Mapping  by 
Numerical  Methods.  National  Bureau  of  Standards  U.  S.  Dept,  of  Commerce  NBS  Tech.  Note 
337,  May  (Reissued  as  ESSATech.  Note  ERL  107  ITS  75  May  1969). 

4  Bowser,  C..  AN/FPS  1 18  Analytical  Model  Specification,  Preliminary  Draft,  Contractor 
General  Electric  Company,  private  communication. 
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3.  MODEL  ANALYSIS 


The  analysis  examined  four  major  areas  of  the  OTH  ionospheric  model:  1)  the  auroral  E 
region  enhancement,  2)  the  midlatitude  trough  region,  3)  the  updating  of  the  OTH 
ionospheric  model  using  observations  from  stations  of  the  AWS  Digital  Ionospheric 
Sounding  Systems  (DISS)  net,  deployed  in  the  vicinity  and  inside  the  coverage  area  of  the 
East  Coast  OTH  radar  operation,  and  4)  the  interpolation  scheme  used  when  data  from  two 
DISS  stations  cU'e  applied  to  update  a  common  area. 


3.1  Auroral  E  Layer 


The  auroral  oval  (Feldstein  and  Starkov  1967)^,  auroral  E  layer,  and  the  F  layer  trough 
boundaries  used  in  the  OTH  ionospheric  model  were  shown  in  Figure  1  of  the  report  by 
Millman  et  al.  (1988).  This  figure  is  reproduced  as  Figure  1  in  this  report  for  convenience. 
Note  that  the  auroral  E  layer  bounded  by  <tifj  (with  the  foE^  maximum  located  at  (|)^| )  and  (|)s 
covers  only  the  equatorward  half  of  the  Feldstein-Starkov  (1967)  oval  defined  by  (|)p  and  (t)E. 
with  as  the  oval  center  line.  For  Q  >  3,  <t>s,  and  (|)£  coincide.  When  Q  <  3,  is  1° 

poleward  of  <t)E.  The  trough  boundaries  are  shown  by  ‘t>Sf  "The  location  of  the 

lowest  trough  density  is  1.5°  equatorward  of  the  northern  trough  boundary  and  has  no 
Kp  dependence.  The  equations  defining  various  boundaries  have  been  added  to  Figure  1. 

The  auroral  foE^  modeled  by  the  OTH  ionospheric  model  is  shown  in  Figure  2  for  0000 
Corrected  Geomagnetic  (CGM)  Local  Time  (LT).  The  CGM  latitude  dependence  of  foEa(MHz)  is 
shown  in  the  figure  for  various  Q  values.  The  horizontal  lines  show  the  width  of  the 
Feldstein-StcU'kov  auroral  oval.  Based  on  the  early  work  of  Wagner®  (1972),  the  OTH 
ionospheric  model  defines  the  auroral  E  layer  contribution  peaking  at  the  equatorial  edge 
of  the  oval.  The  contribution  extends  in  a  triangular  shape,  only  up  to  the  edge  of  the  oval. 
There  is  no  auroral  E  layer  in  the  poleward  region  beyond  the  center  of  the  oval.  The 
auroral  E  layer  affects  the  propagation  situation  of  the  East  Coast  OTH  Radar  System  (ECRS) 
only  for  part  of  the  night.  The  radar  is  >  400  km  due  south  of  the  average  night  time 
southwest  extent  of  the  auroral  E  layer,  which  translates  into  ranges  >  800  km  away  from 
the  radcir  in  the  direction  of  the  northernmost  sectors  of  segment  1 .  Therefore  only  low 
angle  rays  are  likely  to  be  propagated  via  an  auroral  mode.  However,  the  auroral  E  layer  is 
an  important  source  of  ionospheric  spread  clutter  to  the  ECRS. 


®  Feldstein,  Y.I.  and  Starkov,  G.V.  (1967)  Dynamics  of  auroral  belt  and  polar  magnetic 
disturbances.  Planet.  Space  Set  15:209-229. 

®  Wagner,  R.A.  (1972)  Modeling  the  Auroral  E-layer  in  Arctic  Ionospheric  Modeling.  Five 
related  papers  by  Gassman,  G.E.,  Buchau,  J.,  Wagner,  R.A.,  Pike,  C.A.  and  Hurwitz,  M.G.,  Air 
Force  Surveys  in  Geophysics,  No.  241,  AFCRL-72-0305,  AD748796. 
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Figure  1.  Geomagnetic  Orientation  of  the  High  Latitude  Ionosphere. 


WIDTH  OF  AURORAL  OVAL  AT  0000  COLT 


Latitudinal  Dependence  of  foEg  for  Various  Q  Values  From  the  GE  Model. 


Reviewing  the  latitudinal  distribution  of  the  auroral  E  maximum  electron  density. 
Whalen^  (1983),  and  Whalen  et  al  ®  (1985)  have  shown  that  statistically  the  auroral  E  layer 
peaks  at  the  center  of  the  Feldstein-Starkov  oval  and  the  auroral  precipitation  flux  tapers 
off  to  the  edges  of  the  oval,  following  a  Gaussian  distribution.  Using  these  recent  results, 
Eqs.  (5).  (9),  and  (10  )  in  Millman  et  al.  (1988)  have  to  be  replaced  by  the  following  set: 


*t*N  ~  Eq.  (5) 


(1) 


•I’m  ‘t>c 


(2) 


08  ‘t*E  (retained  from  Eq.  (1)  of  Millmcm  et  al.  (1988)) 


(3) 


where  the  L.H.S.  refers  to  the  foE^  boundaries  and  the  R.H.S.  refers  to  the  Starkov^  (1969) 
oval  boundaries.  In  the  OTH  model  the  latitudinal  distribution  of  foE^  is  given  by  Eqs.  (12) 

and  (13): 


foEa=  0.6*(foE3)M+0.4*(foEjM  -<{»m)/(<1>M  "  <t>s)  ^q.  (12)  Millman  et  al.  (1988) 


foEa  =  0.6*(foEa  )m  +0.4*(foEa  )m  ((^mn  "  4>M)/(‘t)N  -  <t>M  )  Eq.  (13)  Millman  et  al.  (1988) 


Using  Whalen's  (1983)  results,  these  two  equations  are  replaced  by  a  single  equation: 


foEa  =  (foE^  exp  (^-3.  *  (^(((t)M  -  <t>)  /  (<t>N  -  4>m  )f 


(4) 


The  results  from  these  modifications  are  shown  in  Figure  3.  Note  that  now  the  peak 
auroral  foE^  is  at  the  center  of  the  oval  and  drops  off  symmetrically  to  the  oval  edges 


7  Whalen,  J.A.  (1983)  A  quantitative  description  of  the  spatial  distribution  and  dynamics 
of  the  energy  flux  in  the  continuous  aurora,  J.  Geophys.  Res.,  88A9:7155-69. 

8  Whalen,  J.A.  O'Neil,  R.R.  and  Picard,  R.H.  (1985)  The  Aurora,  in  Handbook  of  Geophysics 
and  the  Space  Environment,  Air  Force  Geophysics  Laboratory.  AFGL-TR-85-0315,  ADA16700. 

8  Starkov,  G.V.  (1969)  Analytical  representation  of  the  equatorial  boundary  of  the 
auroral  oval  zone.  Geomag.  and  Aeronomy  (Eng.  Ed).  9:614. 
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according  to  the  fourth  root  of  a  Gaussian  distribution.  The  latitudinal  distribution  of 
precipitating  particle  flux  is  a  Gaussian.  The  fourth  root  enters  the  equation  because  the 
relation  between  the  incident  particle  flux  and  the  auroral  foE^  is  given  by 


Flux  «  . 


(5) 


Note  that  in  Figures  2  and  3  the  auroral  foE^  curve  for  Q=0  is  separate  from  the  family  of 
curves  for  other  Q  values.  This  happens  because  in  the  OTH  model  for  Q=0,  the  Feldstein- 
Starkov®  oval  for  Q  is  used,  whereas  for  the  remaining  Q  values  Starkov's  (1969)  empirical 
equation  is  used.  As  a  result  the  oval  moves  4°  equatorward  from  Q=0  to  Q=1  and  at  a  rate 
of  one  degree  in  latitude  for  each  additional  increment  of  AQ=1  from  Q=1  to  Q=9. 

The  auroral  E  layer  from  the  OTH  ionospheric  model  and  the  proposed  changes  are 
shown  in  Figures  4  and  5  for  Q=3.  and  in  Figures  6  and  7  for  Q=6.  In  the  figures  each  block 
presents  the  Corrected  Geomagnetic  Local  Time  (CGMLT)  marked  at  the  left-hand  top  corner. 
The  vertical  lines  show  the  foE^  boundaries.  The  comparison  shows  that  the  present  OTH 
ionospheric  model  uses  only  1/2  of  the  auroral  oval  width,  and  the  auroral  E  distribution 
peaks  at  the  equatorward  boundary  instead  of  at  the  center  as  in  the  improved  Whalen 
(1983)  model.  The  proposed  changes  are  Important  to  the  radar  coverage  when  the  radar 
will  actually  be  close  to  or  under  the  auroral  E  region,  as  it  will  be  in  Alaska. 


3.2  Midlatitude  F  Layer  Trough 

Current  understanding  of  the  behavior  of  the  night  time  trough  suggests  a  region  of  reduced 
electron  density  to  the  south  of  the  equatorward  boundary  of  the  oval,  extending  from  ~  1 
hour  after  sunset  to  the  sunrise,  when  sunlight  covering  the  polar  cap  provides  a  source  of 
ionization,  filling  the  trough  region  by  antisunward  convection.  The  OTH  ionospheric  model 
does  not  properly  reproduce  this. 

(In  his  recent  work,  Whalen  (1987, 1989)  has  shown  the  existence  of  the  trough  wall 
into  the  daytime  ionosphere  along  the  equatorward  boundary  of  the  afternoon  oval.  This 
portion  of  the  trough  has  been  appropriately  named  the  "daytime  trough".  The  extent 
towards  noon  is  strongly  controlled  by  the  magnetic  activity.  Due  to  its  more  northerly 
location  than  the  post  sunset  trough,  it  is  not  relevant  to  the  ECRS  and  has  not  been 
included  in  the  current  considerations.  However,  since  it  may  be  a  source  of  daytime 
ionospheric  clutter  for  the  ECRS,  and  it  may  be  of  substantial  Importance  to  the  ARS 
(Alaska  Radar  System),  we  will  investigate  its  importance  to  the  OTH  radars  at  a  later  date.) 


Whalen,  J.A.  (1987)  The  daytime  F  layer  trough  observed  on  a  macroscopic  scale, 
J.  Geophys.  Res.,  92:2571. 

D  Whalen,  J.A.  (1989)  The  daytime  F  layer  trough  and  its  relation  to  ionospheric  - 
magnetospheric  convection,  J.  Geophys.  Res.,  94:17169. 
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Figure  4.  Latitudinal  Distribution  for  Various  CGT  Hours  for  9=3  From  the  GE  Model. 
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Latitudinal  Distribution  for  Various  CGT  Hours  for  Q=6  From  the  GE  Model. 


The  correct  definition  of  the  trough  is  of  Importance  to  the  ECRS,  since  the  trough  at 
night  can  extend  deep  into  segment  1.  The  trough  plays  an  even  more  important  role  in 
defining  the  propagation  situation  for  the  ARS,  where  it  covers  a  large  part  of  the  ARS 
segment  10  ionospheric  reflection  region  for  18  hours  of  each  day  (moderate  conditions). 

The  ARS  segment  9  will  be  affected  by  the  trough  for  approximately  6  hours  each  night.  The 
equations  defining  the  midlatitude  F  layer  trough  in  the  OTH  ionospheric  model,  which  are 
of  concern  in  the  present  discussion,  are  listed  below. 


AN  =  t{1+cos  [2n  (D+1 1)/365I1 

Eq.  (36)  Millman  et  al.  (1988) 

T  =  Tj  exp  I(Xi-Xi2)/2]  exp  {[-  (t-3)2l/12} 

Eq.  (37)  Millman  et  al.  (1988) 

Xi  =  ((l)-<t»E)/XA 

Eq.  (38)  Millman  et  al.  (1988) 

Xa=  3.7+  1.3  Kp 

Eq.  (39)  Millman  et  al.  (1988) 

tj  =  -0.2 

Eq.  (41)  Millman  et  al.  (1988) 

(h™F2)Mt  =  450-100  V 

Eq.  (44)  Millman  et  al.  (1988) 

R  =  0.25-  l(t-3)2/61  +  2  loge  (-10  tj)  + 

Eq.  (51)  Millman  et  al.  (1988) 

2  loge  (1  +  cos  I27C  (D+1 1)/36511 


The  midlatitude  F  layer  trough  computed  from  the  OTH  ionospheric  model  is  shown  for 
Q=3  and  9=6  in  Figures  8  and  9  respectively  for  Julian  Day  354.  In  the  figures  each  block 
presents  the  dependence  of  the  trough  reduction  factor  of  foF2  with  COM  latitude  for  the 
CGMLT  shown  in  the  left  hand  top  corner.  The  equatorward  boundary  (<t)st).  deepest  point 
((|)^jt)  and  poleward  boundary  (^Nt)  of  the  trough  are  marked  by  triangles.  In  the  OTH  model, 
the  trough,  as  a  night-time  phenomenon,  is  correctly  absent  for  07CX)-1700  CGMLT.  The 
figures  however  show  that: 
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Figure  8.  Latitudinal  Distribution  of  the  Trough  Depletion  Factor  (1+n)  for  Various  CGT 
Hours  for  Q=3  (Kp=l)  From  GE  Model. 
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Figure  9.  Latitudinal  Distribution  of  the  Trough  Depletion  Factor  (1+n)  for  Various  CGT 
Hours  for  Q=6  (Kp=4)  From  GE  Model. 


1)  The  modeled  trough  Is  very  wecik  or  absent  up  to  2300  CGMLT.  This  results  from 
the  vanishing  of  x  (Eq.  (37),  Millman  et  al.  1988)  for  1800  <  t  <  2200  CGMLT. 

Satellite  observations  (Halcrow  eind  Nisbet  1977)  show  that  the  trough  is  formed 
about  11/2  hrs  after  sunset.  Goose  Bay  ionospheric  soundings  also  show  the  onset  of 
the  trough  after  sunset  for  low  magnetic  activity.  At  Goose  Bay  the  onset  time  moves 
to  pre-sunset  time  as  the  magnetic  activity  Increases  (Buchau  1990,  private 
communication).  Recent  eilrcraft  measurements  and  resulting  trough  maps  support 
this.  These  results  are  presented  in  Section  3.3. 


2)  The  modeled  trough  has  a  sharp  knife  edge  with  no  relation  between  the  lowest  value 
of  the  modifying  factor  and  the  latitude  (see  Figure  1).  The  location  of  the  lowest 
point  of  the  trough  can  be  found  by  differentiating  the  latitude  dependence  term 

exp  (-(Xi-Xi2)/2)  in  Eq.  (37)  of  Millman  et  al,  which  yields  the  value  Xj  =  0.5  . 
Substituting  from  Eq.  (39)  in  Millman  et  al,  we  have 

(t.Mt  =  0.5Xa  +  (t)E  =  .5  (3.7+ 1 .3  Kp) 

Thus  contour  of  the  lowest  electron  density,  in  the  trough  lies  equatorward 

of  <t)E.  the  equatorward  boundary  of  the  auroral  oval  (same  as  the  northward 
boundary  of  the  trough)  by  1.85  +  0.65  Kp  degrees  of  latitude.  This  combined  with  the 
definition  of  (jj^t  (=<tiE)  produces  a  northern  vertical  wall  of  the  trough. 

3)  The  lowest  trough  critical  frequency,  is  0.6  foF2.  This  corresponds  to  a  reduction  of 
0.36  in  electron  density  («  foF22'  in  the  trough  region.  Satellite  observations show  a 
reduction  to  0.25  in  the  electron  density  (or  a  factor  of  0.5  in  foF2)  in  the  trough 
region. 

4)  A  look  at  Eq.  (36)  (Millman  et  al.  1988)  shows  that  the  trough  is  absent  in  summer 
(Julian  day  (JD)  170)  and  deepest  in  winter  (JD  354).  When  we  increase  the  Julian  Day 
in  steps  of  10  days  the  first  occurrence  of  the  trough  is  seen  on  JD  220.  Thus  the 
trough  is  absent  from  JD  170  -(220-170)  to  JD  220,  that  is,  the  trough  is  absent  for 
(50*2=)  100  days,  or  more  than  3  months.  This  likely  underestimates  the  occurrence 
of  the  trough  in  summer  and  C£m  be  readily  adjusted.  We  are  currently  analyzing 
Goose  bay  data  to  improve  the  occurrence  statistics. 


12  Halcrow,  B.W.  and  Nisbet,  J.S.  (1977)  A  Model  of  F2  Peak  Electron  Densities  in  the  Main 
Trough  Region  of  the  Ionosphere,  Radio  Science.  12:815-820. 
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5)  Figures  8  and  9  show  a  reduction  of  the  trough  width  at  2300  hours.  To  understand 
this  reduction  in  width,  the  CGMLT  dependence  of  the  trough  width  is  determined 
from  Eqs.  (34)  and  (48)-(50)  (Millman  et  al.  1988),  and  is  shown  in  Figure  10.  For 
1800-2200  CGMLT  the  trough  has  a  fixed  width  of  1.5°  (and  an  insignificant  depth  as 
seen  in  Figures  8  and  9).  From  2300-0600  CGMLT  the  trough  width  shows  a 
systematic  time  dependence,  provided  by  Eq.  (48)  (Millman  et  al.  1988).  Equation  (48) 
(Millman  et  al.  1988)  which  defines  (jigt  and  Eq.  (51)  (Millman  et  al.  1988)  which 
defines  R,  provide  a  time  dependent  width  for  2300  to  0600  CGMLT.  At  2300  CGMLT 
this  width  is  smaller  than  1.5°,  used  for  prior  hours.  Equation  (49)  (Millman  et  al. 
1988)  assigns  the  width  to  the  trough,  but  does  not  control  the  depth  of  the  trough. 

The  resulting  width  of  the  trough  for  1800-2200  CGMLT  shown  in  Figure  10  has  no 
meaning  since  in  this  time  sector  the  depth  factor  is  zero  or  close  to  zero  (see  also 
discussion  in  1)  above). 


To  overcome  the  problems  listed  above  the  following  changes  are  proposed: 

a)  The  absence  of  the  trough  for  100  days  can  be  reduced  to  60  days  by  using  the  square 
root  of  the  term  {l+cos(27t(D+l  1)/3651)  in  place  of  the  linear  term. 

b)  The  trough's  deepest  point  is  aligned  with  ((>^4  replacing  the  (Xj-Xi2)/2  term  by 
-Xj2/2  term. 

c)  As  the  trough  is  known  to  have  a  steeper  foF2  gradient  (Pike  1976)  on  the  poleward 
side  than  on  the  equator  side,  we  define  the  range  of  the  poleward  wall  as  1/5  the 
width  of  the  trough.  Therefore  we  have 


<l*Mt  “  ‘l*Nt  “  ( ^  /5)(<tiNt  "‘I’st  ^ 


(7) 


The  shape  of  the  trough  is  defined  using  the  following  equations: 


t  =  Ti  exp(-Xi2/2  )  expl-(t-3)2/121 


(8) 


where  Xj  is  defined  by  the  equation 
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Xj  =  2((t)  -  <t)Mt  )/^A 


(9) 


and  Xj  is  replaced  by  2Xi  for  <t>  >  resulting  in  a  steeper  foFa  gradient  for  the 
poleward  trough  wall  than  the  equatorward  trough  wall. 


d)  To  generate  the  trough  for  the  1800-0000  CGMLT  period,  the  t-3  term  is  set  to  zero  for 
this  interval.  This  produces  a  trough  with  constant  depth  and  width  for  the  1800- 
0000  CGMLT  interval. 

e)  Since  the  maximum  depth  of  the  trough  is  estimated  to  be  0.5  of  the  non- trough  foF2, 
Tj  is  adjusted  from -1/5  to -1/3 


Xi  =  -l/3. 


(10) 


The  trough  resulting  from  the  proposed  modification  is  shown  in  Figures  11  and  12 
for  Kp=3  and  6  respectively.  The  figures  show  that:  a)  the  trough  is  now  present  at 
night  for  1800-0600  hours  (CGMLT).  b)  the  trough  depth  varies  from  20  percent  of 
unmodified  foF2  at  1800  CGMLT  to  approximately  50  percent  at  0300  CGMLT,  and 
recovers  to  20  percent  at  0600  CGMLT  (sunrise),  and  c)  the  trough  has  poleward  and 
equatorward  walls;  with  the  poleward  wall  steeper  than  the  equatorward  wall, 
supporting  the  observational  evidence. 


Equation  (44)  (Millman  et  al.  1988)  fixes  the  trough  height  at  an  altitude  450  km  at  0300 
CGMLT  with  smooth  lowering  to  350  km  at  the  terminators.  The  trough  altitude  is  estimated 
to  be  50  km  above  the  ambient  non-trough  height.  Therefore  we  propose  the  trough  altitude 
as  the  height  from  the  ITS-78  model  +50  km.  which  ties  the  height  of  the  trough  to  the 
background  model  height. 
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Figure  11.  Modified  Trough  for  Day  354  Using  OTH  Kp  vs  Q  Relation  for  Q=3  (Kp=l). 


Figure  12.  Modified  Trough  for  Day  354  Using  OTH  Kp  vs  Q  Relation 


3.3  Defining  the  foPj  Trough  From  Observations 


Observations  have  shown  that  the  OTH  ionospheric  model  foF2  trough  definition  is 
inaccurate.  The  absence  of  the  foFa  trough  from  1800-0000  CGMLT  is  not  supported  by  the 
observations.  On  the  other  hand,  the  observations  show  that  the  foF2  trough  is  well 
developed  within  1  to  2  hours  after  sunset.  An  example  of  the  foF2  trough  derived  from  the 
data  recorded  during  an  Airborne  Ionospheric  Observatory  (AlO)  flight  on  2/3  March  1989, 
and  from  the  simultaneous  ground  based  observations,  is  shown  in  Figure  13.  The  data  are 
presented  in  the  COM  Local  Time/Latitude  coordinate  system.  The  foF2  contours  show  that 
the  trough  develops  abruptly  at  1930  CGMLT  and  is  present  at  least  up  to  0230  CGMLT,  the 
end  of  the  AIO  flight.  The  cross  sections  of  the  trough  along  CGM  meridians  are  shown  in 
hourly  increments  from  the  onset  of  trough  development  until  0000  CGMLT  (Figure  14).  The 
figure  shows  that  the  trough  is  well  developed  at  1930  CGMLT.  As  time  progresses,  the  trough 
in  general  shows  equatorward  movement  from  1930-0000  CGMLT.  The  trough  wall  on  the 
poleward  side  (RHS)  is  steeper  than  that  on  the  equatorward  side  (LHS).  The  movement  of 
the  equatorward  wall  is  larger  than  that  of  the  poleward  wall.  The  trough  is  deeper  when  it 
is  formed  and  becomes  shallower  with  the  general  decay  of  the  night  time  foF2.  Trough  cross 
sections  for  two  additional  days  are  shown  in  Figures  15  and  16.  These  data  and  other 
available  data  sets  clearly  support  the  suggested  modification  of  the  OTH  ionospheric  model 
F  layer  trough. 


3.4  K-AWS  /9e  Control  of  Trough  and  Auroral  Oval 

The  F  layer  trough  of  the  OTH  ionospheric  model  is  controlled  by  the  90  minute 
geomagnetic  activity  index  K-AWS,  provided  in  90  minute  increments  by  GWC.  The  auroral 
oval  geometry  and  the  auroral  E  and  F  layer  enhancements  are  controlled  by  an  effective 
oval  diameter  index  Qg,  derived  at  GWC  from  satellite  particle  precipitation  data. 

Relationships  between  the  magnetic  Q  index  (for  the  discussion  here  equivalent  to  Qe)>  Kp 
and  K-AWS  have  been  developed  by  AWS  and  Dandekari^  (1982).  These  equations  and  the 
corresponding  graphs  are  shown  in  Figure  17.  The  AWS  set  of  equations  is  incorporated  in 
the  OTH  ionospheric  model,  to  permit  the  Qg  dependent  oval  diameter  (equatorward  oval 
boundary)  to  be  determined  from  K-AWS  in  the  absence  of  the  GWC  Qg  data  that  is,  in  the 
absence  of  the  satellite  data. 


13  Dandekar,  B.S.  The  statistical  relations  among  Q,  Kp  and  the  Global  Weather  Central 
K-indices,  Environmental  Research  Papers,  No.  763,  AFGL-TR-82-0010,  ADAl  18734. 
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2-3  March  1989  From  AIO  Flight  and  Ground 
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Figure  14.  Longitudinal  Cross-Sections  of  2-3  March  1989  foF2  Trough  for  Several  Hours. 


GROUND  BASED 


Figure  16.  Longitudinal  Cross-Sections  of  7-8  March  1989  foF2  Trough  for  Several  Hours. 
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Figure  17.  Kp  vs  Q  Relation  From  AFGWC  Eind  Dandekar  (1982).  Solid  line  -  least  squares  fit. 


The  current  GWC  procedure  for  the  generation  of  from  satellite  particle  precipitation 
data  uses  a  set  of  CGM  local  time  dependent  equations  of  the  precipitation  boundary 
parameterized  in  Kp  (Gussenhoven  and  Hardy  1983;^'^  AFGL  Handbook  of  Geophysics  and  the 
Space  Environment,  Chapter  12,  1985).  When  Qg  is  not  available  from  the  particle 
precipitation  data,  it  is  derived  from  the  ground  based  K-AWS  index  (available  at  90  minute 
Intervals  using  the  Hardy  (private  communication)  algorithm).  The  auroral  boundary 
displayed  at  the  ECRS/WCRS  in  the  E/A  model  maps  and  displays  the  precipitation 
boundary,  determined  from  the  Gussenhoven-Hardy  algorithms  rather  than  the  Starkov 
oval  boundary  used  in  the  OTH  ionospheric  model.  Its  location  is  provided  as  a  set  of 
coordinate  pairs  in  the  AWS  messages.  The  oval  coincides  with  the  precipitation  boundary 
only  in  the  midnight  sector.  In  all  other  sectors  the  oval  boundary  is  to  the  north  of  the 
precipitation  boundary,  farthest  in  the  noon  sector  (~5°).  (However,  it  should  be  noted  that 
the  OTH  ionospheric  model  uses  the  oval  boundary  (Starkov  equation)  as  the  demarcation 
line  between  the  trough  and  the  auroral  E  and  F  layers.  The  AWS  provided-Qg  is  the 
parameter  controlling  the  oval  boundary.)  A  limited  set  of  Qg/K-AWS  data  available  from 
several  (DT&E)  campaigns  is  shown  in  Figure  17  together  with  a  least  squared  deviation  fit 
straight  line.  When  the  data  points  superpose,  the  number  on  the  right  hand  side  lists  the 
number  of  occurrences.  The  data  clusters  and  the  straight  line  fit  very  close  to  that  of  the 
Hardy  algorithm  Indicate  that  most  of  these  Q  data  are  determined  from  the  ground  based  K- 
AWS  observations  using  the  Hardy  algorithm  and  not  from  the  satellite  borne  particle 
precipitation  data.  It  is  clear  that  substantial  differences  exist  between  the  three  sets  of 
equations  relating  the  oval  diameter  index  Qg  and  K-AWS.  Our  observations  also  Indicate  a 
potential  systematic  overestimation  of  the  oval  diameter.  We  are  currently  Investigating 
GWC  procedures  and  the  underlying  assumptions  and  will  provide  an  analysis  of  the 
problem  in  the  necir  future. 

In  Table  2  in  Millman  et  al.  *  the  auroral  latitudes  listed  under  2000  CGM  time  are  too 
low  by  2°.  This  error  has  occurred  because  the  original  Feldstein-Starkov^  curves  for  this 
interval  were  plotted  2°  too  far  south  (drafting  error).  The  interpolation  scheme  used  in  the 
model  propcigates  the  error  proportionally  over  the  1630-2330  CGMLT  sector.  Therefore  a  2° 
correction  has  to  be  added  to  the  latitude  values  of  2000  CGMLT  and  the  corresponding 
corrections  should  be  applied  for  1630  to  2330  CGMLT  intervals  in  Table  2  of  Millman  et 
al.^ 


3.5  Updating  the  OTH  Ionospheric  Model  With  Real  Time  Vertical  Incidence  (VI) 
Site  Data 

The  specification  of  the  ionosphere  by  the  OTH  ionospheric  model  is  improved  by 
updating,  using  real  time  data  from  several  Digisonde  stations  deployed  in  or  close  to  the 
OTH  coverage  area.  The  procedure  involves  an  interpolation  scheme.  As  this  information  is 


Gussenhoven,  M.S.,  Hardy,  D.A.  and  Heinemann,  N.  (1983)  Systematics  of  the  equator- 
ward  diffuse  auroral  boundary,  J.  Geophys.  Res.,  88:5692. 
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not  available  In  the  published  report  (Mlllman  et  al.)^  we  used  the  OTH  supporting  document 
(unpublished  private  communication,  Carl  Bowser,  GE)  for  this  discussion. 

The  Interpolation  scheme  Is  based  on  a  correlation  study  by  Rush  and  Glbbs^^  (1973) 
using  data  from  N-S  and  E-W  (geomagnetic)  chains  of  Ionospheric  stations.  The  weight 
factors  used  in  the  interpolation  scheme  are  the  correlation  coefficients,  which  are  q=l  at 
the  Vertical  Incidence  (VI)  sounder  sites  (separation  0  nmi)  and  which  decrease  to  q=0  as  the 
separation  between  the  VI  sites  Increases  to  a  certain  range,  at  which  the  ionospheric 
variations  at  the  two  sites  become  uncorrelated  (independent ).  For  foF2  and  M(3000)  Rush 
and  Gibbs  (1973)  observed  that  for  a  given  station  these  correlation  distances  are  different 
In  the  E-W  and  N-S  direction  and  they  show  a  very  strong  seasonal  dependence. 

In  the  OTH  ionospheric  model,  the  base  model  is  the  ITS-78  (Barghausen  et  al.^®  1969) 
which  needs  the  sunspot  number,  Julian  day  and  Universal  Time.  To  adjust  the  ITS-78  for 
specification  purposes  AFGWC  determines  an  effective  sunspot  number  SSN^jy  from  the  past 
five  days'  foF2  observations  from  a  network  of  50  VI  ionospheric  stations.  This  SSN^jf  is 
used  to  drive  the  ionospheric  model  providing  a  best  fit  to  the  global  data.  On  the  base  model 
auroral  E  and  F  layer  enhcmcements  (function  of  Q  index)  and  F  layer  trough  (function  of 
Kp)  are  superposed.  The  resulting  model  is  the  AFGWC  polar  ionospheric  model.  A  second 
correction  is  then  applied  that  uses  real  time  data  from  local  sounder  stations  in  the 
vicinity  of  the  coverage  area. 

The  interpolation  scheme  uses  the  model  parameters  predicted  for  a  given  hour.  It 
computes  the  percent  change  between  the  observation  and  the  prediction  of  a  given 
parameter  at  the  ionospheric  station  providing  the  data.  To  the  grid  points  in  the  vicinity  of 
that  site  a  fraction  of  the  observed  change  is  applied.  The  fraction  is  computed  as  the 
correlation  coefficient  (Rush  and  Gibbs),*®  which  is  a  function  of  the  separation  distance 
between  the  VI  site  and  the  grid  point.  In  this  way  all  the  grid  points  in  the  coverage  area  are 
updated.  Such  updated  predictions  based  on  real  time  observations  would  indeed  be  more 
realistic  than  the  simple  'effective  sunspot  SSN^jy '  driven  model  result. 

In  using  this  approach  the  OTH  ionospheric  model  suffers  from  several  errors. 


(1)  Presentation  of  correlation  coefficients.  As  the  weight  factors  are  different  along  the 
E-W  and  N-S  directions  (Rush  and  Gibbs  1973),  the  shape  of  the  weight  factors  as  a 
function  of  distance  between  the  grid  point  and  the  VI  site  will  be  an  ellipse.  The 
curve  determined  from  Eq.  (90-41)  (p.  1878  unpublished  report  by  Carl  Bowser)  is 
shown  in  Figure  18.  Two  items  to  be  noted  in  Figure  18  are:  a)  the  weight  factor  peaks 
at  a  certain  distance  away  from  the  VI  site  (located  at  (0.0))  and  b)  the  curve  is  not 


*®  Rush,  C.M.  and  Gibbs  (1973)  Predicting  the  day-to-day  variability  of  the  mid-latitude 
ionosphere  for  applications  to  HF  propagation  predictions.  Air  Force  Surveys  in  Geophysics 
No.  268.  May,  AFCRL-TR-73-0335,  AD76471 1. 

*®  Barghausen,  A.F.,  Finney,  J.W.,  Proctor,  L.W.  and  Schultz,  L.D.  (1969)  Predicting  long¬ 
term  operational  parameters  of  high  frequency  sky  wave  telecommunication  systems,  ESSA 
Technical  Report  ERL  1 1 0-rTS-78,  Washington,  DC. 
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elliptical.  In  the  computer  program  the  values  greater  than  1  are  set  equal  to  1.  but 
this  does  not  alter  the  shape  of  the  curves  in  Figure  18. 


(2)  Assignment  of  areas  of  updating  to  Vertical  Incidence  (VI)  sites.  The  scheme  uses  at 
most  two  stations  for  updating  a  given  grid  point  and  assigns  somewhat  arbitrary 
cireas  to  the  individual  VI  sites.  The  VI  sites  used/planned  for  the  OTH  update  are 
listed  in  Table  1.  Figure  19a  shows  the  OTH  coverage  area  and  the  locations  of  the  VI 
sites  currently  included  in  the  updating  scheme.  Additional  Digital  Ionospheric 
Sounding  System  (DISS)  stations  of  Importance  to  the  East  Coast  Radar  System 
(ECRS)  model  update  are  now  available  since  the  deployment  of  DISS  at  Sondrestrom 
and  Narsarssauq,  Greenlemd  and  at  Ramey  AFB,  Puerto  Rico.  The  OTH  coverage  area 
is  broken  into  segments  (1-3)  and  sectors  (1-8  in  each  segment)  for  which  the 
Identified  stations  provide  primary  update  (Figure  19a).  Of  these  the  Azores  site  will 
not  be  established,  while  Patrick  AFB  and  Croughton,  UK  have  not  yet  become 
operational.  Goose  Bay.  Argentia.  and  Bermuda  are  inside  the  OTH  coverage  area 
while  Bangor.  Wallops  Island,  and  Patrick  AFB  are  on  the  border  of  the  coverage 
area.  The  closest  grid  point  of  the  radar  coverage  area  is  450  nm  away  from  the  VI 
site  Croughton.  UK.  For  Croughton  the  program  uses  1-1/2  hour  old  data  to  update 
foFa  in  the  coverage,  but  fails  to  transfer  the  VI  site  to  the  appropriate  (22.5°  West  of 
Croughton)  location,  into  the  coverage  area  of  the  radar.  As  it  is.  it  introduces  an 
unwanted  error  due  to  the  diurnal  variation  at  Croughton  over  the  1-1/2  hour  time 
interval.  Figure  19b  shows  the  area  allocation  of  secondary  VI  sites.  The  allocation 
of  the  grid  points  to  the  VI  sites  used  by  the  OTH  model  from  Table  XLIl  p.l882  (Carl 
Bowser  unpublished  report),  for  updating  are  listed  in  Table  2.  The  primary  sites  are 
not  underlined.  The  secondary  sites  are  underlined.  The  grid  point  closest  to  the  VI 
site  is  marked  by  Note  that  the  sites  2.  4.  and  6  are  used  as  secondary  in  their  own 
immediate  vicinity. 


Table  1.  Vertical  Incidence  (VI)  Update  Sites  for  East  Coast  Radar  System  Coverage 


NO 

STATION  NAME 

GEOGRAPHIC 

LAT.  LONG. 

0° 

0 

BANGOR  USA 

44.98 

291.18 

3.9 

1 

THULE.  GREENLAND 

76.50 

291.33 

31.4 

2 

GOOSE  BAY.  CANADA 

53.31 

299.40 

7.9 

3 

CROUGHTON.  UK 

52.00 

358.83 

16.4 

4 

ARGENTIA.  CANADA 

47.29 

306.03 

8.4 

5 

SAN  MIGUAL  ISLAND.  AZORES 

38.00 

335.00 

11.5 

6 

HAMILTON.  BERMUDA 

32.33 

295.33 

3.4 

7 

WALLOPS  ISLAND.  USA 

38.00 

284.50 

1.7 

8 

PATRICK  AFB.  USA 

28.17 

279.33 
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Figure  18.  Weight  Factor 
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Figure  19a.  OTH  Coverage  Area,  Showing  tJie  Locations  of  Vertical  Incidence  Sites  and  Regions  Allotted 
to  the  Respective  VI  sites  Where  That  Site  Will  be  the  Primary  Source  for  Updating  Ionospheric  Parameters. 


Figu: 


Figures  19a  and  19b  and  Table  2  show  that  the  assignment  of  the  areas  to  the 
respective  DISS  stations  is  somewhat  arbitraiy.  While  the  present  eirea  assignments 
prevent  the  stations  in  the  north  from  influencing  areas  in  the  south  (that  is.  trough 
stations  cannot  influence  the  midlatitude  ionosphere),  it  also  prevents  the  smoothing  of 
the  resulting  model  by  the  Joint  use  of  the  available  updating  data  base  during  times 
when  all  stations  are  in  a  very  similar  ionosphere,  namely  during  daytime,  when  all 
stations  are  under  the  influence  of  the  midlatitude  ionosphere. 

Special  problems  arise  with  the  current  scheme,  if  the  updating  data  base  is 
Incomplete.  Non-availability  of  data  from  certain  stations  (for  example,  Argentia, 
station  4)  will  lead  to  unreliable  updates  (fallback  to  five  days’  average)  over  large 
Irregular  areas  with  possible  shcirp  transitions  at  the  borders. 

A  look  at  Table  2  shows  that  only  32  percent  of  the  grid  points  are  updated  by  two  VI 
sites,  while  a  large  fraction  (68  percent)  have  to  rely  on  a  single  update. 


(3)  Update  using  two  Vertical  Incidence  (VI)  sites.  The  two-site  weight  factors  K'j  and  K'2 
listed  in  Eqs.  (90-45,  46)  are  defined  on  pp  1879-80  (Carl  Bowser  unpublished  report). 
These  are: 


K'i  = 


1-q 


>2 

12 


(90-45) 


q  q  *  q 

^  02~  ^  01  *  ^  12 


1-q 


>2 

12 


(90-46). 


The  use  of  an  incorrect  equation  for  the  weight  factor  (refer  to  Figure  18)  yields 
Incorrect  values  for  q'op  q'o2-  and  q'i2-  Also,  when  the  denominator  l-q’^12 
approaches  zero,  these  equations  can  yield  large  positive  and  negative  values  for  K’j. 

The  K'j  values  computed  from  this  GE  model  can  be  positive  and  negative  as 
shown  in  Table  3.  In  Table  3  the  left  hand  column  lists  the  segment  and  the  sector 
(beam)  number  of  a  grid  point.  The  remaining  columns  at  the  top  refer  to  the  grid 
point  number  along  a  ray.  The  first  grid  point  is  200  nmi  from  the  transmitter  site. 
The  consecutive  beams  are  separated  by  7.5°  in  azimuth.  Notice  in  columns  10  and  11 
that  the  first  two  weight  factor  pairs  have  comparable  magnitudes  but  opposite  signs. 
Also  along  segment  1  sector  5,  the  grid  points  5  and  7  have  K'  values  of  0.502,  but  the 
grid  point  number  6  between  the  two  has  the  K'  value  2.120.  Such  discontinuous 
values  result  in  errors  in  the  update. 

To  avoid  an  infinite  value  of  K'  it  is  arbitrarily  set  to  0.5  when  the  denominator 
in  Eqs.  (90-44),  and  (90-45)  is  zero,  in  the  respective  computer  program. 
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Table  3.  Weight  Factor  K  for  Two  Station  Update 


o 

00 

cs 

csi 

in 

o 

0) 

CO 

00 

CO 

d 

d 

f-H 

(N 

CD 

m 

in 

CD 

M 

o 

C'J 

CD 

d 

d 

d 

d 

in 

00 

00 

N 

03 

CO 

03 

03 

d 

d 

d 

00 

CD 

o 

CD 

m 

o 

CO 

d 

d 

d 

o 

2 

D 

'-3 


m 

2 

o 

td 

W 

I 

D 

I 

M 


g 

b: 

8 

1 

O 

2 

P 

K 

O 


If) 

g 

c? 

.s 


o 

CO 

in 

!  00 

Q 

CO 

03 

03  :  : 

i?5 

2) 

03 

03 

q 

•  d 

d 

d 

d 

d 

in  in 


o 


?  ? 


O  -H 

<N  in 
^  in 

d  d 


00 


. 2^^® 

a  o . z^in-HM 


.9 

M 

9 

in 

d 

in 

d 

in 

CO 

o 

in 

(N 

GO 

d 

d 

d 

.  (N 

:  S  S 

•  d  d 


03 

C'J 

o 

§ 

cs 

in 

m 

s 

m 

s 

CO 

QO 

N 

in 

CO 

9 

m 

9 

d 

d 

d 

d 

CO 

CO 

CO  :  : 

.  o 

.  ,  C'J 

<N . 

S  ;  ;  ;  :  : 

. 03 

.  .  .  .  .  ^ 

g  a 

^  CO 

CO 

in 

a 

m 

a 

in 

a 

fS 

,  .  csi 

d . 

. C'J 

1 

9  9 

d 

d 

d 

o 

o 

.  CJ 

(N . 

.  9 

. in 

N 

03 

CO 

cs 

03 

03 

in 

m 

m 

m 

CSI 

in  !  1  ;  !  a 

S  ! 

. 00 

q 

M 

a 

a 

iR 

ih 

2 

•  d 

d . 

. 9 

1 

1 

V 

d 

o 

d 

o 

o 

o  o  o 


o 

d 


(N 

CN 

T^ 

on 

n 

q 

CO 

o 

03 

in 

in 

cs 

CN 

CO 

o 

d 

d 

d 

o 

O 

o 


m  in  S 

m  in  9 

d  d  d  7  7  ?  7 


00 

a 

03 

N 

o 

cs 

CO 

cs 

oo 

in 

00 

CD 

03 

8 

m 

a 

in 

a 

in 

a 

CO 

1 

cs 

c4 

d 

9 

d 

d 

d 

d 

o  o  o  o 


M 

03 

03 

in 

CO 

in 

a 

CO 

N 

CO 

03 

C'J 

CO 

C'J 

a 

d 

oj 

d 

d 

o 

d 

o 

03 

cs 

h- 

CD  03 

00 

CO 

in 

CO 

N  CD 

CO  N 

25 

q 

03 

00 

CO  00 

’’t  a 

o 

dddd<??'?77<?<?77 


S 

O 

? 


in  o 
o  § 
?  O 


fc 

M 

CO 

in 

CO 

N 

00 

M 

CO 

in 

CD 

N 

OO 

N 

CO 

in 

g 

h- 

00 

2 

s 

2 

2 

2 

S 

S 

2 

2 

2 

s 

2 

S 

s 

S 

s 

S 

s 

S 

S 

S 

0 

CLi 

CD 

CD 

CQ 

CQ 

CQ 

CQ 

CQ 

m 

m 

m 

CQ 

CQ 

CD 

CQ 

CQ 

CQ 

CD 

CQ 

CQ 

03 

QQ 

CQ 

CQ 

CQ 

Q 

(N 

(N 

CJ 

IN 

cs 

N 

CD 

CD 

CO 

CD 

CO 

CO 

CO 

CO 

S 

o 

o 

o 

o 

o 

u 

u 

o 

CJ 

o 

U 

o 

O 

U 

O 

o 

U 

O 

a 

O 

g 

g 

g 

g 

o 

w 

cn 

cn 

CO 

cn 

cn 

cn 

cn 

w 

CO 

cn 

OT 

w 

CO 

CO 

cn 

W 

W 

cn 

W 

cn 

cn 

cn 

cn 

§  8 


S 


36 


OTOO- 


These  weight  factors  were  proposed  by  Gautier  eind  Zacharisen^^  (1985)  for 
predicting  variations  at  a  given  station  from  the  observations  at  the  other  station  (or 
at  a  later  time  for  the  same  station)  provided  the  seasonal  and  diurnal  standard 
variations  at  each  hour  at  both  stations  are  known.  Since  these  standard  variations 
are  not  known  for  the  pairs  of  the  OTH  grid  point  and  the  VI  site,  these  equations 
should  not  be  used. 

Corrected  values  for  K'j  were  computed  from  the  above  equations  with  the  proper 
system  of  coordinates  (geographic  in  place  of  geomagnetic)  used  for  computation  of 
distance  of  a  grid  point  from  a  VI  site  and  the  use  of  an  ellipse  (see  next  section) 
instead  of  Eq.  (90-41)  (p.l878  Carl  Bowser  unpublished  report)  for  weight  factors. 

With  this  approach  the  joint  weight  factors  K'j  and  K'ii  for  the  stations  ’i'  and  'ii'  show 
a  slow  reduction  with  increasing  distance  of  the  grid  points  from  the  paired  VI  sites 
and  one  of  the  values  (either  K'j  or  K’„  but  not  both)  may  be  negative,  but  for  the  given 
pair  the  difference  between  the  two  is  positive.  An  explanation  for  this  is:  depending 
on  the  separation  between  the  pair  of  updating  VI  sites  the  redundancy  of  correction 
by  each  site  is  reduced  (negative  K')  by  an  appropriate  fraction.  We  are  proposing  use 
of  all  available  sites  for  updating  instead  of  the  two  station  limit  of  the  present 
approach.  The  K'  factors  need  to  be  computed  for  multiple  station  updates.  Gautier 
and  Zacharisen*^  have  presented  an  elaborate  method  for  computing  K'  factors  for  a 
multiple  station  update.  A  rather  simple  method  for  replacing  K'j  is  proposed  in  the 

following  section. 

(4)  Separation  between  Vertical  Incidence  (VI)  site  and  the  grid  point.  To  determine  the 
separation  between  the  VI  site  and  grid  point  the  OTH  model  uses  the  geomagnetic 
coordinate  system,  which  results  in  incorrect  separation  distances,  leading  to 
incorrect  weight  factors  (from  Eq.  (90-41)  pi 878,  Carl  Bowser  unpublished  report). 
This  is  due  to  the  fact  that  the  geomagnetic  coordinate  system  is  distorted  with 
respect  to  the  geographic  coordinate  system. 


4.  PROPOSED  CORRECTIONS/MODIFICATIONS  OF  OTH  MODEL 


To  correct  errors  in  procedure  and  shortcomings  in  the  updating  process,  the  following 
steps  are  suggested: 


1.  Correct  the  presentation  of  the  correlation  function  (ellipse). 


Gautier,  T.N.  and  Zacharisen,  D.H.  (1965)  Use  of  Space  and  Time  Correlation  in  Short¬ 
term  Ionospheric  Predictions,  First  Annual  IEEE  Communications  Convention.  Boulder, 
Colorado,  June  7-9,  pp.  671-676. 
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2.  Allow  all  available  DISS  stations  and  the  radar’s  own  VI  sounder  to  participate  in 
the  updating  process  at  all  the  grid  points  in  the  OTH  coverage  area  ionosphere 
according  to  their  distances  to  the  respective  updating  location  and  the  phenomenon. 

3.  Modify  the  two/multiple  station  updating  procedure  with  a  simple  algorithm. 

4.  To  prevent  sub-oval  stations  from  updating  the  auroral  and  trough  F-layers  and  vice 
versa,  we  propose  to  use  the  Qe  Starkov  equatorward  oval  boundary  and  the  K-AWS 
trough  boundary  to  determine  the  northward  extent  of  influence  of  sub-trough 
stations,  and  the  southward  influence  of  the  oval  stations. 

5.  The  trough  will  be  updated  only  by  stations  like  Goose  Bay  and/or  Argentia,  which 
have  actually  moved  into  the  trough  region  as  defined  by  the  location  of  the 
equatorward  boundary  of  the  oval  and  of  the  equatorward  boundary  of  the  trough. 
The  trough  is  bounded  by  the  1800  (through  0000)  and  0600  CGMLT  meridians. 


Specifically.  Goose  Bay  and/or  Argentia  move  from  the  midlatitude  ionosphere  in  the 
daytime,  into  the  trough  ionosphere  after  sunset,  and  into  the  oval  at  a  time  depending  on 
Qg.  We  propose  here  to  let  Goose  Bay  and  Argentia  update  the  midlatitude  ionosphere  In  the 
daytime  according  to  the  weight  of  the  correlation  function.  As  the  Qg/K-AWS  defined 
oval/and  trough  move  into  the  coverage  area,  the  trough/oval  area  rvlll  not  be  updated  by 
the  DISS  stations,  since  at  this  early  time  no  DISS  will  be  in  or  even  close  to  these  regions. 
Goose  Bay  and  Argentia  will  only  affect  the  trough  or  oval  ionosphere  if  the  specified  trough 
or  oval  boundary  is  south  of  the  respective  station.  The  influence  of  the  stations  on  the 
updating  process  will  then  be  limited  to  the  respective  (trough  or  oval)  regime. 


Presentation  of  correlation  coefficients:  The  Eq.  (90-41)  p.l878  (Carl  Bowser  unpublished 
report)  is  replaced  by  the  equation  of  an  ellipse  using  the  following  procedure.  The  equation 
of  an  ellipse  is  given  by 


(11) 
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where  a  and  b  are  the  major  (E-W)  and  minor  (N-S)  axes  of  the  ellipse. 

Figure  20  illustrates  the  approach  for  computing  the  weight  factor  (or  the  correlation 
distance).  The  figure  shows  a  large  ellipse  computed  from  Eq.  (11).  This  is  the  boundary  at 
which  the  weight  factor  'q'  is  zero.  The  updating  VI  site  is  at  the  center  of  this  ellipse.  The 
grid  point  'P'  to  be  updated  is  at  a  distance  r  and  at  an  azimuth  0  measured  from  the  X-axis, 
which  Is  the  major  axis  of  the  ellipse.  The  line  from  the  center  passing  through  P  cuts  the 
ellipse  at  the  distcmce  'R'.  At  the  center  of  the  ellipse  (that  is,  at  the  VI  site)  'q'  is  1  .  If  r  is 
greater  than  R,  the  point  is  outside  the  boundary  of  the  correlation  ellipse,  therefore  q  is  set 
to  zero  (q=0).  The  value  'q'  at  the  point  P  is  given  by  the  equation 


q  =  l =  l-r*- 

R 


cos^  0  sin^  0 


(12) 


In  computing  a  distance  and  azimuthal  location  between  two  points  on  the  earth  the 
amgle  is  measured  with  respect  to  north.  For  the  ellipse  (Figure  20),  it  is  measured  from  the 
east  with  respect  to  the  major  axis.  Therefore  90°  should  be  subtracted  from  the  azimuthal 
angle  between  the  pair  of  the  stations  . 

The  axes  of  the  ellipse  'a'  emd  'b'  show  seasonal  dependence  for  both  foF2  and  M(3000) 
parameters  (Rush  and  Gibbs  1973).  The  quantities  ’a’  and  'b'  are  constant  for  a  given  season 
for  a  given  parameter  (foF2,  MUF). 

Rush  and  Gibbs  (1973)  determined  the  weight  factors  for  the  geomagnetic  (N-S  and  E-W) 
meridian  reference  frame.  For  this  purpose  the  axes  of  the  ellipse  referred  to  in  Eq.  (12)  have 
to  be  rotated  so  as  to  align  the  minor  (b)  axis  along  the  geomagnetic  N-S  meridian.  These 
angles  are  computed  for  each  station  and  listed  in  the  last  column  of  Table  1 .  Taking  into 
account  a)  the  azimuthal  difference  between  geographic  and  geomagnetic  systems  of 
coordinates,  and  b)  the  eizimuthal  cmgle  for  the  ellipse  with  respect  to  east  vs  the  azimuthal 
angle  between  stations  with  respect  to  north,  the  angle  9  for  Eq.  (12)  is  computed  by  the 
equation 


9=9^  +92  -90° 


(12a) 


where  9i  is  the  angle  from  Table  1  for  the  VI  site  and  02  ‘S  azimuthal  angle  between  the  VI 
site  and  the  grid  point. 


Updating  with  VI  site  data:  To  overcome  the  shortcomings  of  the  current  updating  scheme 
(see  Section  3.5,  item  2  above)  we  propose  to  use  all  the  sites  for  updating  the  model  at  a 
given  grid  point  out  to  the  maximum  range  of  their  respective  correlation  distances  (R=0) 
with  restrictions  in  items  2  and  4.  Therefore  we  replace  Eq.  (90-52)  by  the  equation 
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(13) 


X  ^ - 

LQi 


where  Xj  Is  the  percent  difference  between  the  observation  and  the  prediction  in  the 
ionospheric  model  parameter  at  the  i*^^  VI  site,  Qj  is  the  distance  dependent  weight  factor 
defined  in  Eq.  (12)  and  the  summation  is  done  for  all  the  VI  sites  from  which  the  real  time 
data  are  available  for  the  time  of  updating  the  model.  Note  that  the  correlation 
coefficient/weight  factor  'q,'  is  the  level  of  reliability/confidence  in  the  covariation,  but 
does  not  provide  the  magnitudes  of  the  respective  variations. 

As  discussed  above,  the  updating  scheme  is  proposed  for  all  DISS  and  ECRS  VI  sounders 
in  daytime,  with  limits  after  sunset,  to  prevent  cross-talk  from  trough/ oval  stations  into 
the  midlatitude  ionosphere  and  vice  versa.  The  area  of  influence  limits  are  as  previously 
discussed,  the  Qg/K-AWS  driven  trough  and  oval  boundaries. 

The  seasonally  dependent  constants  'a'  and  'b'  in  Eq.  (12)  are  computed  from  the  data  in 
Table  XLl  pp.  1876-1877  (Carl  Bowser  unpublished  report).  A  straight  line  is  fitted  to  each 
data  set  with  the  condition  the  line  passes  through  (0,1.0),  so  that  the  correlation  is  1.0  at 
the  VI  site  where  the  separation  distance  is  zero.  The  equation  for  this  straight  line  is  given 
by 


y  =  mx+ 1 . 


(14) 


Using  the  least  squeu'ed  deviation  technique,  we  get  the  slope  m  from  the  equation 


m  = 


Xyn-n 


(15) 


where  the  summation  is  done  for  all  the  ’n’  non-zero  (weight)  data  points. 

The  results  for  the  straight  line  fit  to  the  data  in  Table  XLl  (pp  1876-1877,  Carl  Bowser 
unpublished  report)  are  presented  in  Table  4. 
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Table  4.  Empirical  Straight  Line  Fit  of  Weight  Factors  to  Data  in  Table  XLl 
(p.  1876-77  Carl  Bowser  unpublished  report) 


PARAMETER 

SEASON 

NOV-FEB 

MAR-APR 

SEPT-OCT 

MAY-AUG 

DIRECTION 

foF2 

E-W 

-0.271 

-0.122 

-0.143 

N-S 

-0.350 

-0.203 

-0.210 

M(3000) 

E-W 

-0.308 

-0.318 

-0.460 

N-S 

-0.580 

-0.351 

-0.713 

The  correlation  coefficients  and  the  straight  line  fit  are  shown  in  Figures  21-23  for  the 
parameter  foFg  and  in  Figures  24-26  for  M(3000).  In  each  figure  we  have  two  sets  of  weight 
factors:  a)  along  the  N-S  direction,  and  b)  along  the  E-W  direction.  Note  that  each  straight 
line  passes  through  the  point  (0,1).  The  cutoff  along  the  distance  axis  is  the  value  for  the 
constant  'a'  for  the  E-W  direction  and  'b*  for  the  N-S  direction.  The  'a'  and  'b'  are  the  major 
£md  minor  axes  of  the  ellipse  in  Eq.  (12). 

The  equal  weight  contours  for  these  data  are  computed  using  Eq.  (12).  These  are  shown  in 
Figures  27-29  for  foF2  and  in  Figures  30-32  for  M(3000).  Note  that  the  ellipse  is  largest  for 
foF2  for  the  Apr-May,  and  Sept-Oct  periods  and  is  smallest  for  M(3000)  for  the  period  May- 
Aug.  One  may  compare  the  shape  of  curves  in  Figures  23-28  with  that  shown  in  Figure  18 
determined  from  Eq.  (90-41)  on  p.l878  (Carl  Bowser,  unpublished  report ). 

In  the  computer  one  can  now  use  three  sets  (for  three  seasons)  of  'a'  and  'b  values  listed 
in  Table  5  for  each  of  the  parameters  foFa  and  M(3000)  in  place  of  the  data  in  Table  XLl  on 
pp.  1876-77  (Carl  Bowser,  unpublished  report). 


Table  5.  Major  (a)  and  Minor  (b)  Axes  (km)  of  Correlation  Ellipse 


PARAMETER 

SEASON 

a 

b 

foF2 

NOV-FEB 

3690 

2850 

MAR-APR 

SEPT-OCT 

8180 

4920 

MAY-AUG 

7010 

4770 

M(3000) 

NOV-FEB 

3250 

1720 

MAR-APR 

SEPT-OCT 

3150 

2850 

MAY-AUG 

2170 

1400 

42 


CORRELATION  COEFFICIENT 


DISTANCE  (km) 


Figure  21.  Dependence  of  the  Correlation  Coefficient  (Weight  Factor)  on  Separation 
Between  VI  Sites  Along  N-S  and  E-W  Directions  for  foFa  for  November-Februaiy 
Period. 
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CORRELATION  COEFFICIENT 


DISTANCE  (km) 


Figure  22.  Dependence  of  the  Correlation  Coefficient  (Weight  Factor)  on  Separation 
Between  VI  Sites  Along  N-S  and  E-W  Directions  for  foFa  for  March-April  and 
September-October  Period. 
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CORRELATION  COEFFICIENT 


Figure  23.  Dependence  of  the  Correlation  Coefficient  (Weight  Factor)  on  Separation 
Between  VI  Sites  Along  N-S  and  E-W  Directions  for  foF2  for  May-August  Period. 
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CORRELATION  COEFFICIENT 


Figure  24.  Dependence  of  the  Correlation  Coefficient  (Weight  Factor)  on  Separation 
Between  VI  Sites  Along  N-S  and  E-W  Directions  for  M(3000)  for  November-February 
Period. 
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CORRELATION  COEFFICIENT 


Figure  25.  Dependence  of  the  Correlation  Coefficient  (Weight  Factor)  on  Separation 
Between  VI  Sites  Along  N-S  and  E-W  Directions  for  M{3000)  for  March-April  and 
September-October  Period. 
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CORRELATION  COEFFICIENT 


DISTANCE(km) 


Figure  26.  Dependence  of  the  Correlation  Coefficient  (Weight  Factor)  on  Separation 
Between  VI  Sites  Along  N-S  and  E-W  Directions  for  M(3000)  for  May-August  Period. 
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Figure  27.  Weight  Factor  Contours  for  foF2  for  November-February  Period. 
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Figure  28.  Weight  Factor  Contours  for  foF2  for  March-April  and  September-October  Period. 
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Figure  29.  Weight  Factor  Contours  for  foF2  for  May-August  Period. 
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Figure  30.  Weight  Factor  Contours  for  M(3(XX))  for  November-February  Period. 
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Figure  31.  Weight  Factor  Contours  for  M(3000)  for  March-April  and  September-October  Period. 
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Figure  32.  Weight  Factor  Contours  for  MOOOO)  for  May-August  Period. 


The  distance  x  and  y  of  the  grid  point  from  the  VI  sites  is  computed  using  the  geographic 
system  of  coordinates  (not  geomaignetic  coordinates  as  done  in  the  GE  approach)  and  the 
ellipse  is  aligned  along  the  magnetic  meridian  using  the  values  of  the  angle  '0'  between  the 
geomagnetic  and  geographic  meridiem  at  the  VI  site  listed  in  Table  1 . 


5.  CONCLUSIONS 


This  preliminary  report  summarizes  the  results  of  two  tests  of  the  OTH  model  installed 
at  the  PL  facility.  Several  inaccuracies  in  the  concept  and  formulation  of  the  auroral  E 
layer.  F  layer  trough  and  the  updating  scheme  using  real  time  Vl-site  data  have  been 
summarized.  A  way  to  correct  the  respective  inaccuracies  has  been  discussed.  In  the  near 
future  we  plan  to  complete  these  modifications  and  test  the  'improved'  model  with 
additional  data  collected  during  the  OTH  radar  DT&E  campaign  conducted  in  October  1989. 


6.  MISCELLANEOUS  CORRECTIONS  AND  COMMENTS 


As  the  GE  document  (Carl  Bowser  unpublished  report)  is  an  indispensable  source  for  the 
procedure  for  coordinate  registration  and  the  ionospheric  modeling  used  at  the  OTH  radars, 
additional  typographical  and  other  errors  and  comments  are  listed  below.  For  convenience 
and  simplification  these  are  listed  in  sequential  order. 

pn.l863  B.  C.  and  D:  The  equations  use  a  term  600  in  the  denominator.  This  term  should  be 
the  denominator  only  for  the  negative  term  in  the  bracket.  The  computer  program  of  the 
corresponding  part  is  written  properly  with  no  error. 

p.l863  E  line  5:  should  read  mult=W_ccc_Vll-GP*mult;  on  line  3  from  the  bottom  2V1 
should  read  V12. 

p.l863  E.F:  There  are  two  F2  layer  adjustment  schemes;  one  listed  here  and  the  other  on 
pp.  1878-80.  The  scheme  listed  here  is  proper  and  better  than  the  one  on  pp.  1878-80. 
n.l865:  in  Eq.  (90-30  the  first  negative  (-)  sign  on  the  R.H.S.  should  be  replaced  by  a 
positive  (+)  sign. 

p.1870:  in  Eq.  (90-23)  the  positive  sign  (+)  in  the  [)  bracket  should  be  replaced  by  a 
negative  (-)  sign. 

n.l873:  in  Eq.  (90-35)  [Be-4AeCe11/2  should  read  -IBe2-4AeCe1^''2. 
n.l876:  col  3.  line  8.  Table  for  E-W  direction.  .09  should  read  .07. 

P.1876:  bottom  line  columns  4  and  5  should  read  .14  and  .17  respectively. 
p.l877:  column  3.  lines  6.7.8.  Table  of  E-W  direction  should  read  .18.  .05.  and  .02 
respectively. 

P.1877:  last  column,  no.  5.  Table  for  N-S  direction  line  4  should  read  .15. 
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p.1880:  Eq.  (90-52)  The  denominator  2  is  needed  if  both  Kj  and  K2  are  non-zero.  If  either  is 
zero  the  real  value  of  Xp  would  be  divided  by  2. 

P.1884:  lines  7-8,  one  needs  vertical  height  estimates  from  a  vertical  incidence  sounder  in 
the  vicinity  where  the  reflections  of  the  radar  and/or  target  signal  takes  place,  not  in  the 
vicinity  of  the  radar  or  the  target. 

p.l961:  lines  22-23.  This  condition  of  removing  FI  layer  when  foFj  >  foF2  results  in 
discarding  the  low  altitude  Fj  layer  ionization.  The  Fj  layer  is  allowed  to  build  up  to  the 
time  it  becomes  equal  to  or  stronger  than  foF2.  In  effect  we  are  throwing  away  the  stronger 

low  altitude  ionization. 

p.l966  In  Eq.  (190-41)  should  read 

P.1967:  In  Eqs.  (190-49)  and  (190-50)  the  equations  should  have  a  negative  (-)  sign  following 
the  equal  to  (=)  sign. 

p.  1967:  On  line  15  the  equation  numbers  (190-4)  and  (190-6)  quoted  refer  to  the  E  layer  and 
not  to  the  F2  layer  under  discussion. 

P.1968:  On  line  10  the  Eq.  (190-4)  refers  to  the  E  layer  and  not  to  the  F2  layer  under 
discussion. 
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